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Abstract 

A series of tensile relaxation tests is performed on isotactic polypropylene in the 
sub-yield and post-yield regions at room temperature. Constitutive equations are 
derived for the time-dependent response of a semicrystalline polymer at isothermal 
loading with small strains. Adjustable parameters in the stress-strain relations 
are found by fitting experimental data. It is demonstrated that the growth of the 
longitudinal strain results in an increase in the relaxation rate in a small interval of 
strains in the sub-yield domain. When the strain exceeds some critical value which 
is substantially less than the apparent yield strain, the relaxation process becomes 
str ain-indep endent . 



* Corresponding author 



1 Introduction 



This paper is concerned with the nonhnear viscoelastic behavior of semicrystalhne poly- 
mers at isothermal uniaxial deformation with small strains. The study concentrates on 
the time-dependent response of isotactic polypropylene at room temperature, but we be- 
lieve that the constitutive equations to be derived may be applied to other semicrystalhne 
polymers [polyamides, polyethylene, polytetrafluorethylene, poly (ethylene terephthalate) , 
poly(butylene terephthalate), poly(ethylene-2,6-naphthalene dicarboxilate), etc.] at tem- 
peratures above the glass transition point for the amorphous phase. The viscoelastic 
and viscoplastic behavior of polypropylene (PP) has been a focus of attention in the 
past decade, which may be explained by numerous applications of this material in in- 
dustry (oriented films for packaging, reinforcing fibres, nonwoven fabrics, polyethylene- 
polypropylene copolymers, blends of PP with thermoplastic elastomers, etc.). 

The nonlinear viscoelastic response of polypropylene was studied by Ward and Wolfe 
see also 0, and Smart and Williams three decades ago, and, more recently, by 
Ariyama |Q, ^, ^, 0, Wortmann and Schulz ^, Dutta and Edward |]rU|, and Read and 



Tomlins |]rT|. In the past couple of years, the linear viscoelastic behavior of PP fibres 



was analyzed by Andreassen [jT2| and Lopez- Manchado and Arroyo [T^ , and the dynamic 
viscoelastic response of blends of isotactic polypropylene with polyethylene and styrene- 



butadiene-styrene triblock copolymer was studied by Souza and Demarquette [|1J] and 



Gallego Ferrer et al. |T^, respectively. It was found that the loss tangent of isotactic 
polypropylene demonstrates two pronounced maxima being plotted versus temperature. 
The first maximum (/5-transition in the interval between T = —20 and T = 10 °C) 
is associated with the glass transition in the most mobile part of the amorphous phase, 
whereas the other maximum (a-transition in the interval between T = 50 and T = 80 °C) 
is attributed to the glass transition in the remaining part of the amorphous phase. This 
conclusion is confirmed by DSC (differential scanning calorimetry) traces for quenched 
PP that reveal an endoterm at T = 70 °C which can be ascribed to thermal activation of 
amorphous regions with restricted mobility under heating [0. 

Scanning electron microscopy and atomic force microscopy evidence that isotactic PP 
is a semicrystalhne polymer with a complicated morphology. The crystalline phase con- 
tains brittle monoclinic a spherulites (the characteristic size of 100 fim) consisting of 
crystalline lamellae (10 to 20 nm in thickness) directed in radial and tangential direc- 
tions [|r^ , and ductile hexagonal /3 spheruhtes consisting of parallel-stacked lamellae . 
The amorphous phase is located between spherulites and inside the spherulites between 
lamellae. It consists of (i) relatively mobile chains between spherulites and between radial 
lamellae and (ii) severely restricted chains in regions bounded by radial and tangential 
lamellae. Mechanical loading results in inter-lamellar separation, rotation and twist of 
lamellae, fine and coarse slip of lamellar blocks, and their fragmentation. The latter leads 



to reorganization of blocks and strain-induced smectic-monoclinic transitions |16|. Strain 



ing of polypropylene causes chain slip through the crystals, breakage and reformation of 
tie chains, and activation of restricted amorphous regions driven by lamellar disintegra- 
tion. In the post-yield region, these transformations of microstructure lead to the onset 



of voids between lamellae, break-up of crystalls, and formation of fibrills [19 



It is hard to believe that these changes in the internal structure of PP may be ade- 
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quately described by a constitutive model with a relatively small number of adjustable 
parameters. To develop stress-strain relations for polypropylene, we apply the method 



of "homogenization of micro-structure" pO|. According to this approach, an individual 



(equivalent) phase is taken into account whose deformation captures essential features of 
the response of a semicrystalline polymer with a sophisticated morphology. 

An amorphous phase is chosen as the equivalent phase because of the following reasons: 

1. The viscoelastic response of semicrystalline polymers is conventionally associated 
with rearrangement of chains in amorphous regions [|l^ . 



According to the Takayanagi-Nitta concept [^, 0, sliding of tie chains along and 
their detachment from lamellae play the key role in the time-dependent response of 
semicrystalline polymers. 

The viscoplastic flow in semicrystalline polymers is assumed to be "initiated in the 
amorphous phase before transitioning into the crystalline phase" p3[ . 



Conventional models for polyethylene |^ , polypropylene |25| and poly(ethylene 
terephthalate) p6|, ^ treat these polymers as networks of macromolecules. 



Above the glass transition temperature for the mobile amorphous phase, polypropylene 
is thought of as a network of chains connected by junctions. Deformation of a specimen 
induces slip of junctions with respect to the bulk material. Sliding of junctions reflects 
slippage of tie molecules along lamellae and fine slip of lamellar blocks which are associated 
with the viscoplastic behavior of a semicrystalline polymer. 

The viscoelastic response of PP is described in terms of the concent of transient net- 
works [^, and it is modelled as separation of active chains from their junctions 
and attachment of dangling chains to temporary nodes. Unlike previous studies, the net- 
work of polymeric chains is assumed to be strongly inhomogeneous, and it is treated as 
an ensemble of meso-regions (MR) with various potential energies for the detachment of 
active strands. Two types of MRs are distinguished: (i) active domains where strands 
separate from junctions as they are thermally agitated (these MRs model a mobile part 
of the amorphous phase), and (ii) passive domains where detachment of chains from junc- 
tions is prevented (these MRs reflect a part of the amorphous phase whose mobility is 
severely restricted by radial and tangential lamellae). Detachment of active chains from 
the network is treated as a thermally-activated process whose rate is given by the Eyring 



formula |32| with a strain-dependent attempt rate. Deformation of a specimen results in 

(i) an increase in the concentration of active MRs (which is ascribed to a partial release of 
the amorphous phase in passive meso-domains driven by fragmentation of lamellae) and 

(ii) a mechanically-induced acceleration of separation of strands from their junctions in 
active MRs. 

An increase in the rate of relaxation of amorphous glassy polymers at straining is 
conventionally assumed to occur in the entire region of nonlinear viscoelasticity (from the 



strain e of at least 0.005 to the yield strain ey), see It is ascribed to the mechanically- 



induced growth of "free volume" between macromolecules which results in an increase in 



their mobility p3, R3, pq, 37|. Despite wide applications of the free volume concept in 



nonlinear viscoelasticity of solid polymers, experimental verification of this approach was 
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limited to relatively small strains (substantially below the yield point). The objectives 
of this study are (i) to report experimental data in tensile relaxation tests on isotactic 
polypropylene in the range of strains between e = 0.005 (which corresponds to the end of 
the region of linear viscoelasticity) and e = 0.16 (at higher strains necking of specimens 
is observed at the chosen cross- head speed), (ii) to derive constitutive equations for the 
time-dependent behavior of semicrystalline polymers at isothermal uniaxial stretching 
with small strains, (iii) to find adjustable parameters in the stress-strain relations by 
fitting observations, and (iv) to assess the interval of strains where the rate of relaxation 
monotonically increases with strain intensity. 

The exposition is organized as follows. Section 2 is concerned with the description of 
the experimental procedure. Kinetic equations for sliding of junctions and detachment 
of active strands and attachment of dangling strands in active MRs are developed in 
Section 3. The strain energy density of a semicrystalline polymer is determined is Section 
4. In Section 5 stress-strain relations are derived for uniaxial deformation of a specimen 
by using the laws of thermodynamics. These constitutive equations are employed to fit 
experimental data in tensile relaxation tests in Section 6. Some concluding remarks are 
formulated in Section 7. 

2 Experimental procedure 

A series of uniaxial tensile relaxation tests on polypropylene specimens was performed 
at room temperature. Isotactic polypropylene (Novolen llOOL) was supplied by Targor 
(BASF). ASTM dumbbell specimens with length 14.8 mm, width 10 mm and height 
3.8 mm were injection molded and used without any thermal pre-treatment. Our DSC 
measurements (the sample mass 16.43 mg, the heating rate 10 K/min) demonstrated 
the specific enthalpy of melting Aifm = 95.5 J/g, which corresponds to the degree of 
crystallinity = 0.46 (with reference to the enthalpy of fusion for a fully crystalline 
polypropylene is assumed to be 209 J/g). 

Mechanical experiments were performed with a testing machine Instron-5568 equipped 
with electro-mechanical sensors for the control of longitudinal strains in the active zone of 
samples (the distance between clips was about 50 mm). The tensile force was measured 
by the standard loading cell. The engineering stress cr was determined as the ratio of the 
axial force to the cross-sectional area of the specimens in the stress-free state (38 mm^). 

16 relaxation tests were carried out at longitudinal strains in the range from ei = 0.005 
to ei6 = 0.16, which correspond to the domain of nonlinear viscoelasticity, sub-yield and 
post-yield regions for isotactic polypropylene (the yield strain, ey, is estimated by the 
supplier as 0.13). 

Any relaxation test was performed on a new sample. No necking of specimens was 
observed in experiments. In the kth relaxation test {k = 1,...,16), a specimen was 
loaded with the constant cross- head speed 5 mm/min (which roughly corresponded to 
the strain rate eo = 0.05 min^^) up to the longitudinal strain e^, which was preserved 
constant during the relaxation time = 20 min. The engineering stresses, a, at the 
beginnings of relaxation tests are plotted in Figure 1 together with the stress-strain curve 
for the specimen strained up to eie = 0.16. The figure demonstrates fair repeatability of 
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experimental data. 

The longitudinal stress, cr, is plotted versus the logarithm (log = log^g) of time t (the 
initial instant t = corresponds to the beginning of the relaxation process) in Figures 2 
to 4. These figures show that the stress, a, monotonically increases with strain, e, up to 
eo = 0.08 and remains practically constant when the strain exceeds the critical strain e^. 

3 A micro- mechanical model 

A semicrystalline polymer is treated as a temporary network of chains bridged by junc- 
tions. The network is modelled as an ensemble of meso-regions with various strengths of 
interaction between macromolecules. Two types of meso-domains are distinguished: pas- 
sive and active. In passive MRs restricted by radial and tangential lamellae, inter-chain 
interaction prevents detachment of chains from junctions, which implies that all nodes 
in these domains are thought of as permanent. In active MRs, active strands (whose 
ends are connected to contiguous junctions) separate from temporary junctions at ran- 
dom times when these strands are thermally agitated. An active chain whose end slips 
from a junction is transformed into a dangling chain. A dangling chain returns into the 
active state when its free end captures a nearby junction at a random instant. 

Denote by X the average number of active strands per unit mass of a polymer. Let 
Xa be the number of strands merged with the network in active MRs, and Xp the number 
of strands connected to the network in passive MRs. Under stretching some crystalline 
lamellae (that restrict mobility of chains in passive MRs) break, which results in the 
growth of the number of strands to be rearranged. As a consequence, the number of 
strands in active MRs increases and the number of strands in passive meso-domains 
decreases. This implies that the quantities X^, and Xp become functions of the current 
strain, e, which obey the conservation law 

X,(e)+Xp(6)=X (1) 

Separation of active strands from the network and reformation of dangling strands in 
active MRs are thought of as thermally activated processes. It is assumed that detachment 
of active strands from their junctions is governed by the Eyring equation [Q, where 
different meso-domains are characterized by different activation energies, lD, for separation 
of active strands from the network. 

According to the theory of thermally-activated processes, the rate of separation of 
active strands in a MR with potential energy Co in the stress-free state is given by 

where /cb is Boltzmann's constant, T is the absolute temperature, and the pre-factor Fa 
is independent of energy u and temperature T. Introducing the dimensionless potential 
energy 
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where Tq is some reference temperature, and disregarding the effects of small increments 
of temperature, AT — T — Tq, on the rate of detachment, F, we arrive at the formula 



r = raexp(-u;). (2) 

It is assumed that Eq. (2) is satisfied for an arbitrary loading process, provided that the 
attempt rate. Fa, is a function of the current strain, 

ra = ra(e). 

The distribution of active MRs with various potential energies is described by the proba- 
bility density p{u) that equals the ratio of the number, Ng^{e,u!), of active meso-domains 
with energy lu to the total number of active MRs, 

N,{e,u;)^X,{e)p{u;). (3) 

We suppose that the distribution function for potential energies of active MRs, p(to'), is 
strain-independent . 

The ensemble of active meso-domains is described by the function na(i,T, a;) that 
equals the number of active strands at time t (per unit mass) belonging to active MRs 
with potential energy u that have last been rearranged before instant r G [O.t]. In 
particular, na(0, 0, lu) is the number (per unit mass) of active strands in active MRs with 
potential energy a; in a stress-free medium, 

na(0,0,u;) = 7Va(0,u;), (4) 

and na{t, t, uj) is the number (per unit mass) of active strands in active MRs with potential 
energy cu in the deformed medium at time t (the initial time t — corresponds to the 
instant when external loads are applied to a specimen), 

n^{t,t,uj) = N^{e{t),uj). (5) 

The amount 



-{t,T,u;) 



dr 

t=T 



dr 

equals the number (per unit mass) of dangling strands in active MRs with potential energy 
cu that merge with the network within the interval [r, r -|- dr] , and the quantity 

— (t,T, u;)dT 

is the number of these strands that have not detached from temporary junctions during 
the interval [r, t] . The number (per unit mass) of strands in active MRs that separate 
(for the first time) from the network within the interval [t,t + dt] reads 

^''^t,0,u;)dt, 



dt 
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whereas the number (per unit mass) of strands in active MRs that merged with the 
network during the interval [r, r + dr] and, afterwards, separate from the network within 
the interval [t,t + dt] is given by 



'dtdr 



{t, r, uj)dtdT. 



The rate of detachment, F, equals the ratio of the number of active strands that separate 
from the network per unit time to the current number of active strands. Applying this 
definition to active strands that merged with the network during the interval [r, r + dr] 
and separate from temporary junctions within the interval [i, t -\- dt], we find that 



(6) 



Changes in the function na(t, 0,a;) are governed by two processes at the micro-level: (i) 
detachment of active strands from temporary nodes, and (ii) transition of passive meso- 
domains into the active state under loading. The kinetic equation for this function reads 



dn N dt 

-Q^{t,0,u;) = -r{e{t),u;)n,{t,Q,u;) + -^{e{t),u;)-{t). 

The solution of Eq. (7) with initial condition (4) is given by 



r* ON.. 



- ^ r(€(s),uj)ds 
Jo 



+ 



Jo 



dt 



j T{t{s),uj)ds 



dr. 



It follows from Eq. (6) that 

dn. 



dr 



(t, r, a;) = (^(r, a;) exp — / T[t{s),uj)ds 



where 



t=T 



To determine the function ip{t,u!), we use the identity 

ft Qfi 

n^{t,t,u}) ^n^{t,0,u) + -^{t,r,u})dT. 
Equations (5) and (11) imply that 



n^{t,0,uj) 



dr 



it,r,uj)dr = NMt),uj). 



Differentiating Eq. (12) with respect to time and using Eq. (10), we obtain 



^(t,.) + ^(^,0,a;) + §^it,r,.)dr = ^(e(t), a;)$(t). 



de 



dt 



(7) 



(8) 

(9) 
(10) 

(11) 

(12) 



This equality together with Eqs. (6), (7) and (11) results in 

rt Qji 



(p{t,u;) = T{e{t),uj) 



dr 



{t, T, Uj)dT 



na(t,0,a;) + 

= r(e(t),cu)na(t,t,cu). (13) 
Substituting expression (13) into Eq. (9) and using Eq. (5), we arrive at the formula 



dr 



— j T{e{s),uj)ds 



(14) 



The kinetics of rearrangement of strands in active MRs is described by Eqs. (2), (3), (8) 
and (14). These relations are determined by (i) the distribution function p{uj) for active 
MRs with various potential energies cu, (ii) the function ra(e) that characterizes the effect 
of strains on the attempt rate, and (iii) the function 



«a(e) 



X ' 



(15) 



that reflects mechanically-induced activation of passive MRs. 

Rearrangement of strands in active MRs reflects the viscoelastic response of a semicrys- 
talline polymer. The viscoplastic behavior is associated with the mechanically-induced 
slippage of junctions with respect to their positions in the bulk material. 

Denote by e-a{t) the average strain induced by sliding of junctions between macro- 
molecules (the subscript index "u" means that eu(t) is associated with the residual strain 
in a specimen which is suddenly unloaded at instant t). The elastic strain (that reflects 
elongation of active strands in a network) is denoted by edt). The functions ee{t) and 
e-a{t) are connected with the macro-strain e{t) by the conventional formula 



€{t)^€^{t) + e^{t). 



(16) 



We adopt the first order kinetics for slippage of junctions with respect to the bulk material, 
which implies that the increment of the viscoplastic strain, de^, induced by the growth 
of the macro-strain, e, by an increment, de, is proportional to the absolute value of the 
stress (7, 

c?e„ , . / de\ ^^^^ 



—— = B\a\ sign((T— I, 



de 



dt^' 



where the pre-factor i? is a non-negative function of stress, strain and the strain rate. 



B 



de 



The last multiplier in Eq. (17) determines the direction of the viscoplastic flow of junc- 
tions. It is convenient to present Eq. (17) in the form 



^{t) = B (a (t) , 6 (t) , ^) I a (t) I sign [a(t ) | (t) 



de\ , 



de 



dt' 



dt 



de , , 



e,(0) = 0. 



(18) 



It follows from Eq. (18) that the rate of sliding vanishes when the strain is constant. 
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4 The strain energy density 



Any strand is modelled as a linear elastic solid with the mechanical energy 

where fi is the average rigidity per strand and e is the strain from the stress-free state to 
the deformed state. 

For strands belonging to passive meso-domains, the strain e coincides with eg. Mul- 
tiplying the strain energy per strand by the number of strands in passive MRs, we find 
the mechanical energy of meso-domains where rearrangement of chains is prevented by 
surrounding lamellae, 



(19) 



With reference to the conventional theory of temporary networks ||3T|, we assume that 
stresses in dangling strands totally relax before these strands merge with the network. 
This implies that the reference (stress-free) state of a strand that merges with the network 
at time r coincides with the deformed state of the network at that instant. For active 
strands that have not been rearranged until time t, the strain e{t) coincides with ee(t), 
whereas for active strands that have last been merged with the network at time r G [0, t], 
the strain e{t, r) is given by 

e(t,r) = ec(t) - ec(r). 

Summing the mechanical energies of active strands belonging to active MRs with various 
potential energies, u, that were rearranged at various instants, r G [0,t], we find the 
mechanical energy of active meso-domains. 



* dria 



(t,r,u;) ee(t) - ec(r) 



dr 



(20) 



The mechanical energy per unit mass of a polymer reads 

W{t) = W,{t) + Wp{t). 

Substituting expressions (19) and (20) into this equality and using Eq. (16), we arrive at 
the formula 



W{t) 



2 
+ 



Xp(e(t))(e(t)-6,(t))'(t) + 



du 



na(t,0,cu)(e(t) -eu(t) 



97 



(t,r,.;)((e(t) -e,(t)) - (e(r) -e.(r)))V 



(21) 



Differentiation of Eq. (21) with respect to time results in 



dW 
~dt 



(t) = ^iA{t) 



+ 2Mo(t), 



(22) 
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where 



A{t) = Xp(e(t))[e(i)-e,(t) 

rt 



+ duj^n^{t,0,u}) e{t) - e^{t) 



de dt 

rt 



+ l^^{t,T,u;)[{e{t)-e^{t)) - (6(r) - 6„(r))]'dr|. 



(23) 



Bearing in mind Eqs. (5) and (11), we transform the first equahty in Eq. (23) as follows: 



A{t) 



/ 1 r 

Xp{e{t))+ NMt),u;)du e{t)-e,{t) 
Jo L 

* dn 



roo 



duj 



lo dr 



it,T,u;) e(T) - eu(r) 



dr. 



This formula together with Eqs. (1) and (3) implies that 

r n On r 

A{t) = X[e{t) - e^{t)\ - duj -^{t,T,u;)[e{T) - e^ir) 

Substitution of expressions (6) and (7) into the second equahty in Eq. (23) yields 



dr. 



Ao{t) 



dXr 



de 



'0 de 



(e(i ) , co) du] ^ (t) [e{t) -e^{t)f- A, (t) , 



where 



+ 



r{e{t),uj)duj^n^{t, 0, Lj) [e{t) - e^{t) 
1^ ^{t,r,u;)[{e{t) - e^{t)) - (e(r) - e^{r)) 



dry 



It follows from Eqs. (1), (3) and (25) that 

Ao{t) = -Ai(t). 
This equality together with Eq. (22) results in 



dW 
"df 



where 



, de,, 



A,{t)^2A{t)^{t). 



(24) 



(25) 



(26) 



(27) 
(28) 
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5 Constitutive equations 



For uniaxial loading with small strains at the reference temperature Tq, the Clausius- 
Duhem inequality reads ^ 



^ dQ , , dW , , 1 , .de , . 



(29) 



where p is mass density, and Q is the entropy production per unit mass. Substitution of 
expression (27) into Eq. (29) implies that 



^ 'ait) - ppAit)] + \ [Mt) + Mt) 



P 



> 0. 



(30) 



Because Eq. (30) is to be fulfilled for an arbitrary program of straining, e = e(t), the 
expression in the first square brackets vanishes. This assertion together with Eq. (24) 
results in the stress-strain relation 



a{t) = ppA{t) 



E 



e{t) - e,(t) 



Or 



{t,T,Uj) 



e r 



eu(r) 



dr 



(31) 



where 



E = p^X 

is an analog of the Young modulus. It follows from Eqs. (18), (28) and (31) that 



(32) 



According to Eqs. (26) and (32), the functions Ai{t) and A2{t) are non-negative for an 
arbitrary program of loading, which, together with Eq. (31), implies that the Clausius- 
Duhem inequahty (30) is satisfied. 

Substitution of Eqs. (3), (14) and (15) into Eq. (31) results in the constitutive 
equation 



E 



X 



e(t) -e,(t) -KMt)) / Pi^)duJ 



T{e(t),uj) exp 



— / T{e{s),iu)ds 



e r 



dr 



(33) 



Given functions pi^uj), ra(e) and Ka(e), the time-dependent response of a semicrystalline 
polymer at isothermal uniaxial loading with small strains is determined by Eqs. (2), (18) 
and (33). For a standard relaxation test with the longitudinal strain e°. 



e(t) 



0, t<0, 
e°, t > 0, 



these equations imply that 



|l-/€a(e°)^"p(^) l-exp(-ra(eO)exp(-cu)t) du^ 
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where e° is the strain induced by shding of junctions. Introducing the notation 

Ci(e°) = E{e' - e°), C,{e') = E{e' - el)K^{e'), (34) 
we present this equahty as follows: 



a{t, e°) = Ci(e°) - ^^(e") J^^ p{uj) [l - exp(-r,(e°) exp(-.;)t 
To fit experimental data, we adopt the random energy model with 



du. (35) 



p{uj) = po exp 



2S2 



(uj > 0), p{uj) = (to < 0), (36) 



where fl and S are adjustable parameters, and the pre-factor pq is determined by the 
condition 

p{uj)du} = 1. (37) 







Given a strain e°, Eqs. (35) and (36) are determined by 5 material constants: 

1. the average potential energy for rearrangement of strands Q, 

2. the standard deviation of distribution of potential energies E, 

3. the attempt rate for separation of strands in active MRs from their junctions Fa, 

4. the coefficients Ci and C2. 

Our purpose now is to find these parameters by fitting experimental data. 

6 Fitting of observations 

We begin with matching the relaxation curve measured at the strain eg = 0.025. Because 
the rate of rearrangement. Fa, and the average activation energy, Q, are mutually depen- 
dent [Eqs. (35) and (36) imply that the growth of fl results in an increase in Fa], we set 
Fa = 1 s and approximate the relaxation curve by using 4 experimental constants: Q, S, 
Ci and C2. To find these quantities, we fix the intervals [0, f^max] and [0, Smax], where the 
"best-fit" parameters fl and S are assumed to be located, and divide these intervals into 
J subintervals by the points Qi = iA^ and T,j = jAs {i, j = 1, . . . , J) with = f^max/ J, 
As = Smax/</- For any pair, the integral in Eq. (35) is evaluated numerically 

(by Simpson's method with 200 points and the step A^ = 0.1). The pre-factor po is 
determined by Eq. (37). The coefficients Ci = Ci{i,j) and C2 = C2{i,j) that minimize 
the function ^ 

Jihj) = 5l[^cxp(im) - (^num{tm) , (38) 

where the sum is calculated over all experimental points tm, are found by the least- 
squares method. The stress ciexp in Eq. (38) is measured in the relaxation test, whereas 
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the quantity CTnum is given by Eq. (35). The "best- fit" parameters fl and S minimize the 
function J on the set 

{^i, (hj = l,■■■,^)}■ 
After determining the "best- fit" constants, and E^, this procedure is repeated for the 
new intervals and [Ej_i, E^+i] to ensure good accuracy of fitting. Figure 2 

demonstrates fair agreement between the experimental data and the results of numerical 
simulation with Q = 5.03 and S = 3.05. 

To approximate relaxation curves at other strains, e^, we fix the constants Q and 
E found by matching experimental data at eg and fit every relaxation curve by using 3 
adjustable parameters: Fa, Ci and C2. To find these quantities, an algorithm is applied 
similar to that employed to match the relaxation curve at eg = 0.025. We fix the interval 
[0,Fjnax]; where the "best-fit" attempt rate Fa is supposed to be located, and divide this 
interval into J subintervals by the points Fj = iAp {i = 1, . . . , J) with Ap = Fmax/-^- 
For any Fj, we calculate the integral in Eq. (35) numerically and find the coefficients 
Ci = Ci{i) and C2 = C2{i) minimizing the function (38) by the least-squares method. 
The "best-fit" attempt rate is determined from the condition of minimum for the function 
J' on the set |Fj (z = 1, . . . , J)|. After finding the "best-fit" value, Fj, this procedure 
is repeated for the new interval [Fj„i,Fj+i] to ensure an acceptable accuracy of fitting. 
Figures 2 to 4 show good agreement between the observations and the results of numerical 
analysis. 

For any longitudinal strain e^, the attempt rate, Fa(efc), is determined by matching an 
appropriate relaxation curve. The fraction of active MRs, Ka,{ek), is found from Eq. (34), 

/ N C'2(efe) 

These quantities are plotted versus strain e in Figures 5 and 6. The experimental data 
are approximated by the phenomenological equations 

log Fa = 7o + Jie, Ka,^ko + he, (39) 

where the coefficients 7^ and ki are found by the least-squares method. Figures 5 and 6 
reveal two different regimes of the nonlinear viscoelastic behavior of polypropylene. In 
the region of relatively small strains (less than = 0.02), the quantities log Fa and Ka 
increase linearly with strain, e. In the region of strains exceeding the threshold value, e^,, 
the attempt rate grows rather weakly, and the fraction of active MRs remains constant. 
This observation implies that the free- volume theory may be applied to the description 
of the nonlinear viscoelastic response of polypropylene at relatively small strains only. 
When strains exceed the threshold value, e*, relaxation of longitudinal stresses becomes 
strain-independent. 

This conclusion results in two questions of interest: 

1. What is a micro-mechanism for transition from the nonlinear (strain-dependent) to 
the linear (strain-independent) viscoelastic behavior of isotactic polypropylene at 
the threshold strain, e*. 
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2. Whether other semicrystalhne polymers demonstrate the nonhnear time-dependent 
response in a hmited range of strains only, i.e., whether relaxation of stresses 
in semicrystalhne polymers (unlike amorphous glassy polymers) becomes strain- 
independent at threshold strains, e*, that are noticeably lower than the yield strain, 
ey. 

These issues will be the subject of a subsequent study. 

7 Concluding remarks 

Constitutive equations have been derived for the time-dependent behavior of semicrys- 
talhne polymers at isothermal loading with small strains. To develop stress-strain rela- 
tions, a version of the mean-field approach is employed: a complicated micro-structure of 
isotactic polypropylene is replaced by an equivalent transient network of macromolecules 
bridged by junctions (physical cross-links, entanglements and crystalline lamellae). The 
network is assumed to be strongly inhomogeneous, and it is thought of as an ensemble 
of meso-regions with various potential energies for separation of strands from temporary 
nodes. 

The viscoelastic response of a semicrystalhne polymer is ascribed to the processes of 
detachment and reformation of chains in active meso-domains. Rearrangement of active 
strands is modelled as a thermo-mechanically activated process whose rate obeys the 
Eyring formula. 

The viscoplastic response is described by slippage of junctions with respect to their po- 
sitions in the bulk material. The rate of sliding of junctions is assumed to be proportional 
to the macro-stress in a specimen. 

The mechanical energy is determined as the sum of the strain energies of active strands. 
Constitutive equations arc derived by using the laws of thermodynamics. These relations 
are applied to study relaxation of longitudinal stresses at uniaxial tension of specimens. 

A series of relaxation tests have been performed on isotactic polypropylene at room 
temperature. Adjustable parameters in the stress-strain equations are found by fitting 
observations. The following conclusions are drawn from the analysis of experimental data: 

1. Despite a complicated morphology of isotactic polypropylene, its time-dependent 
response is rheologically simple in the sense that the relaxation spectrum (which is 
determined by distribution function, p(a;), for potential energies of active MRs) is 
not affected by mechanical factors. 

2. Two regimes of the nonlinear viscoelastic behavior of polypropylene are observed. 
At relatively small strains (less than the threshold value = 0.02), the attempt 
rate. Fa, and the fraction of active MRs, k.^, increase with strain, e, in agreement 
with the free- volume concept. When the strain, e, exceeds its threshold value, e*, 
the attempt rate. Fa, grows with strain rather weakly, while the fraction of active 
MRs, Kg,, remains constant. 
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Figure 1: The longitudinal stress a MPa versus strain e. Unfilled circles: experimental data 
in a tensile test with the cross-head speed 5 mm/min. Filled circles: experimental data at the 
beginnings of relaxation tests 
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Figure 2: The longitudinal stress a MPa versus time t s in a tensile relaxation test at a strain 
e. Circles: experimental data. Solid lines: results of numerical simulation. Curve 1: e = 0.005; 
curve 2: e = 0.0075; curve 3: e = 0.010; curve 4: e = 0.0125; curve 5: e = 0.015; curve 6: 
e = 0.0175; curve 7: e = 0.020; curve 8: e = 0.025 
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Figure 3: The longitudinal stress a MPa versus time t s in a tensile relaxation test at a strain 
e. Circles: experimental data. Solid lines: results of numerical simulation. Curve 1: e = 0.03; 
curve 2: e = 0.04; curve 3: e = 0.06; curve 4: e = 0.08 



26.0 



a 



16.0 




1.0 logt 3.0 



Figure 4: The longitudinal stress a MPa versus time i s in a tensile relaxation test at a strain 
e. Circles: experimental data. Solid lines: results of numerical simulation. Curve 1: e = 0.10; 
curve 2: e = 0.12; curve 3: e = 0.14; curve 4: e = 0.16 
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Fi gurc 5: The attempt rate s versus strain e in tensile relaxation tests. Circles: treatment 
of observations. Solid lines: approximation of the experimental data by Eq. (39). Curve 1: 
^0 = -0.66, 71 = 34.73; curve 2: 70 = 0.05, 71 = 0.17 
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Figure 6: The fraction of active MRs Ka versus strain e in tensile relaxation tests. Circles: 
treatment of observations. Solid lines: approximation of the experimental data by Eq. (39). 
Curve 1: ko = 0.40, ki = 6.28; curve 2: ko = 0.49, ki = 0.0 
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